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ABSTRACT 


The  sensitivity  of  planetary  waves  to  various  initial 
conditions  and  thermal  forcing  is  examined  using  a  linear 
global  primitive  eguation  spectral  model. 

Initial  conditions  are  obtained  by  switching  on  a 
analytic  heat  source  and  then  integrating  the  model  equa¬ 
tions  out  to  30  days.  An  averaging  procedure  is  used  to 
eliminate  any  transient  modes  which  remained  after  integra¬ 
tion  so  that  the  initial  conditions  represent  steady  state 
solutions.  Additional  integrations  are  performed  in  which 
•errors*  are  introduced  into  the  forcing  and  initial 
conditions. 

Results  of  the  study  are  examined  using  polar  phase  vs 
amplitude  plots  (harmonic  dials)  of  various  spherical 
harmonics.  Results  indicate  that  planetary  waves  are  not 
sensitive  to  errors  in  the  forcing  or  initial  conditions. 
However,  this  lack  of  sensitivity  is  most  likely  due  to  the 
simplified  initial  conditions  used  in  the  model  rather  than 
to  the  inherent  sensitivity  of  planetary  waves. 
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I.  INTRODUCTION 


Atmospheric  predictability  studies  (Lorenz,  1969)  indi¬ 
cate  each  scale  of  motion  has  a  limit  to  it’s  predict¬ 
ability.  Small  scale  motions  are  theoretically  predictable 
to  an  hour,  synoptic  scales  to  a  few  days  and  planetary 
scales  to  a  few  weeks.  However  a  number  of  studies,  Lambert 
and  Merilees  (1978),  Baumhefner  and  Downey  (1978),  and  Morse 
(1983),  have  shown  that  synoptic  scale  motions,  rather  than 
planetary  scale  motions,  are  the  most  accurately  forecast. 
Fig.  1.1  is  taken  from  Daley  et  al  (1981).  It  gives  the 
500mb  geopotential  forecast  error  from  the  6-level  National 
Center  for  Atmospheric  Research  global  circulation  model  as 
a  function  of  wavenumber.  The  error  for  each  wavenumber  has 
been  normalized  by  the  natural  atmospheric  variance  for  that 
wavenumber.  Note  that  as  previously  stated  the  synoptic 
scales  (wavenumbers  4-8)  are  in  fact  the  most  accurately 
forecast.  This  result  does  not  seem  to  be  dependent  on  the 
type  cf  model  used.  Baumhefner  and  Downey  examined  a  number 
of  different  models  and  found  similar  results.  The  error  in 
forecasting  planetary  scale  waves  might  not  seem  to  be  that 
important  given  that  the  error  in  the  synoptic  scales  is 
smaller  and  that  these  are  the  scales  of  motion  which 
produce  most  of  the  day  to  day  weather  changes.  The  impor¬ 
tance  of  an  accurate  planetary  scale  comes  to  light  when  one 
starts  to  consider  medium  to  long  range  forecasts  (up  to  10 
days).  Since  the  planetary  waves  often  act  to  steer  the 
smaller  synoptic  scale  disturbances  an  improved  planetary 
scale  forecast  would  presumably  lead  to  an  improved  forecast 
on  the  synoptic  scale.  In  addition,  the  planetary  waves 
contain  a  major  portion  of  the  eddy  kinetic  energy  in  the 
atmosphere  and  again  a  better  long  wave  forecast  would  lead 
to  a  better  synoptic  scale  forecast. 
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There  are  at  least  two  possible  reasons  why  planetary 
waves  are  not  forecast  as  well  as  theory  suggests.  One  is 


Figure  1. 1  Forecast  Errors  for  1,  2 
and  3  Days  of  Integration. 

that  the  model  dynamics  for  planetary  waves  are  inadeyuate. 
Given  that  planetary  waves  are  quasi-stationary  in  nature 
and  are  forced  in  part  by  differential  heating  and  orography 
it  would  seem  reasonable  to  suspect  that  inadequacies  exist 
in  the  dynamics  of  the  model’s  forcing.  Another  possibility 
is  that  current  initialization  procedures  are  not  adequate 
for  planetary  scales. 

It  is  the  purpose  of  this  study  to  examine  the  nature  of 
planetary  waves  and  to  determine  the  sensitivity  of 
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planetary  waves  to  initial  conditions  and  forcing.  Chapter 
two  of  this  thesis  describes  the  model  used  in  this  study. 
The  third  chapter  outlines  the  nature  of  the  experiment.  In 
the  fourth  chapter  the  results  of  the  experiment  are 
presented  and  in  the  final  section  these  results  are 
discussed  and  some  conclusions  are  drawn. 


II.  MODEL  DESCRIPTION 


The  model  used  in  this  study  is  a  baroclinic  spectral 
transform  model  which  was  developed  by  fiosmond  (1977). 
Similar  formulations  have  been  carried  out  by  Hoskins  and 
Simmons  (1975)  and  by  Bourke  (  1974).  The  model  is  config¬ 
ured  to  include  friction  and/or  diabatic  heating.  The 
specifics  of  how  friction  and  diabatic  heating  are  included 
in  the  model  will  be  discussed  in  the  next  section.  The 
basic  equations  of  the  model,  in  sigma  coordinates,  are  as 
follows: 


a _c 

3  t 
3D 

at 

ia 

3  t 

39 
3 1 

3 

3a 


-V*  (  ;  +  f  )V-k*^x(  RT^q+d-~)  +ic-Vx? 
k-Vx(  s  +  f  )V-V*  (RTVq+a^-)-V2(4>+f  )+V*F 


-D-V  • 


-v*ye+a|^-+S^ 

3a  PCp 


RT 

a 


where : 

£  -  vorticity 
D  -  divergence 
T  -  temperature 
0  -  potential  temperature 
tt  -  surface  pressure 
V  -  horizontal  velocity  vector 
<t>  -  geopotential  height 
R  -  gas  constant 

Cp  -  specific  heat  at  constant  pressure 
f  -  Coriolis  parameter 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 
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o 


vertical  coordinate  (o  =  p/ir) 

o  -  vertical  velocity  (<5  =  ^2.  \ 

dt 

q  -  In 

p  -  p 
-  R/Cp 

F  -  frictional  force 

The  continuity  equation  (Eg.  2.3)  may  be  rewritten  by 
integrating  with  respect  to  sigma  and  applying  the  boundary 
conditions  a(0)  =  a (1)  =0.  Thus  the  integral  of  Eg.  2.3  maybe 
written  as 


I'l  - 

at  " 


-D+G 


(2.6) 


where  (  )  =  /J(  )da  and  G=  V  *Vg. 

The  vertical  velocity,  <j,  may  be  obtained  diagnostically 
by  substituting  Eg.  2.6  into  Eg.  2.3  and  integrating  in  the 
vertical  to  obtain 


6  =  (D+G)a-/J(D+G)da  '2,7> 

which  uses  a  (0)  =0 

The  first  law  of  thermodynamics  Eg.  2.4  can  be  written 

IS-  =  -v4T-oK6^-(To_<)+<T(^+V-^q)+S-  <2*8> 

o  t  ao  at  cp 

In  order  to  apply  semi-implicit  differencing  it  is  necessary 
to  divide  temperature  as  follows: 


T  =  T*(o)+T’(o,  X,  t) 

where  T*  represents  an  appropriately  averaged  temperature. 

The  basic  equations  can  be  conveniently  written  in 
spherical  coordinates  by  defining  the  following  operator: 
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ct(a,b)  5  -Jy 

l-x*  3X  3x 


(2. 10) 


Osing  Egs.  2.9  and  2.10  the  basic  equations  can  be  written 
as  follows: 


||  =  -a(A,B) 


(2.  11) 


'N  T-N 

3^-  =  -a(A,B)-?2(E+4>+RT*q) 


(2.  12) 


=  -aan”,VT’)+Dr'-aKa|-(Ta-<)+KT(G-G-D) 


(2.  13) 


It  =  “(D+G) 


(2.  14) 


“I?1  -*1 


(2.  15) 


where 

B  =  <?*f  jv-dg-  ♦Si  rk 

°*dl rM 


E  =  (U2+Vz )/2(l-x2) 


U  =  u(cos4')/r 


V  =  v(cosY)/r 


x  =  sin^ 


f  =  latitude 


X  =  longitude 


Eys.  2.10  -  2.15  are  the  basic  equations  used  in  the  model. 
These  equations  are  represented  spectrally  in  the  horizontal 
and  finite  differenced  in  the  vertical. 

A.  VEBTICAL  STRUCTOBE 

The  vertical  structure  of  the  model  follows  the  develop¬ 
ment  given  by  Arakawa  and  Suarez  (1983).  The  variables  are 
staggered  in  a  so  that  D#  0,  7  and  T  are  carried  at  the 

A 

mid-point  of  each  layer  where  a  =  and  a  is  carried  at  the 
top  and  bottom  of  each  layer  where  o  =a^  .  The  vertical 
structure  is  illustrated  in  Fig.  2.  1.  The  finite  difference 
form  of  Eqs.  2.10  -  2.15  are 


3t  =  ak(A,B) 


=  ak(B,-A)-V2(<j)k+RT*q  +  Ek) 


=  - ( G+D) 


=  -ak(UT' ,VT' )T'Dk+  T* ( Gk~G-D) - 

fe^k+lBk(^+^k+1-Tk)+  Ak(Tk-  ^Jk-l)+Q/CP 

k  k+ 1 

*Ltf  VCPT"im<£  -D 


LM'Pt 


a ,  .  =  a,  x1(G+D)-S  (G.+D.)Aa. 

k+1  k+1  3=1  3  3  3 

where 

l'k  =  <ok*i,r)'' 


D  -  /  V  .  , 1  +  K  1  +  K  \  /  / 

l+<)(  k  +  l°k  )/(ak+l"ak) 


(2.  16) 
(2.  17) 
(2.  18) 
(2.19) 

(2.  20) 
(2.21) 
(2.22) 


15 


A  -  <C  +  f>V(4ktlCVkn-V*ak<Vk-Vk_1))/2Ao, 


tRT'k(l-x2  >/rM§§  -(/T^!/r)F( 


L"  • 

K 

B  = 

(I,f)Vliku|uM-V,i 

S  ’ 

■ 

G  = 

-CCpRT^/r)|a  +( /l-x2 /r)F 
V  k  *  ^71-  /  IT 

•  V 

*v 

« 

D  = 

<4 

• 

<4 

7T 

H 


i 

> 

& 


A,  =  (P,  -P,  )/(P.  t,-P,  ) 
k  k  k  k+1  k 


Bk  -  (P 


k+l-Pk)/(Pk+l-P^)=l-Ak 


Egs.  2.16  -  2.22  can  be  written  in  matrix  form,  such 
that  the  terms  on  the  right  hand  sides  contain  all  the  terms 
which  are  to  be  evaluated  explicitly  and  the  left  hand  side 
contains  those  terms  which  are  to  be  evaluated  implicitly. 
Eg.  2.20  can  be  combined  with  an  integrated  finite  differ¬ 
ence  fora  of  Eg.  2.  15  to  obtain 


<p  =  C  T+<J> 


(2.23) 


where  C  is  a  sguare  matrix  and  the  other  quantities  are 
column  vectors.  The  finite  difference  form  of  the  surface 
pressure  tendency  equation (Eg.  2.14)  is 


||  =  I(Gk+Dk)Ack 

which  can  be  written  in  matrix  form  as 


=  -NT(G+D) 

o  t 


(2.  24) 


(2.25) 
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T'/j-'r  '  s’** 

cl*.  V  VuJ  -  t  -  tU A- W  . 


T 

where  N  is  the  transpose  of  a  constant  column  vector. 
Similarly  Eg.  2.22  can  be  written 


o  =  Z( G+D) 

The  next  to  last  term  in  Eg.  2.  19  is 


(2. 26) 


6  B  (fJS—T 

ktl  k  pxn 


ktl‘Tk)tlikAk-l(Tk"  F^k-l' 


(2.27) 


For  the  purpose  of  se  ni-implicit  formulation,  the  tempera¬ 
ture  is  separated  according  to  Eg.  2.9.  The  mean  part  of 
that  term  can  be  written  as 


ak*iBklFfrpk*i-TkHaAk-i<Tk-  =  S(G+D) 

K+l  k-1 

Egs.  2.17,  2.19,  2.23,  and  2.25  may  now  be  written: 


(2.  28) 


($» +RT*q)=a(B,-A)-^2E  =  i?D 
|^+NTD  =  -nt2 


{2. 29) 
(2.  30) 
(2.31) 


(2.32) 


T 

where  <2  =  M  ♦  <T*H  and  4»=  <f>  -  <f>g  and  KD  and  KT  represent 

terms  which  have  not  been  explicitly  separated  out. 

The  semi-implicit  time  differencing  is  achieved  by  eval¬ 
uating  the  terms  on  the  left  hard  sides  of  Egs.  2.29,  2.30 

and  2.31  implicitly.  The  remaining  terms  and  Eg.  32  are 
evaluated  explicitly  using  leapfrog  differencing.  The 


difference  eguations  can  now  be  written 


(2.  33) 
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.  VVr 


-K  I  rw  'JH  ’  ~  W  v.  •»  «•  . 


?n  +  l+At^n  +  l  =  Vl-Atg5n-l-2At(^T)n 
qn  +  l+At^  =  ^n-l_At^T^n-l"2At^T^n 


(2.34) 


(2. 35) 


Now  the  following  eguation  for  D  can  be  found  by  substi¬ 
tuting  Eg.  2.34  and  2.35  into  2.33: 

B  S+,  =  B.S  , +2 At ( K- )  -2AtV2 ( CT  ,+RT*q  ,+C(iL)  (2.36) 

-  n+i  +  n-1  D  n  =  n-1  —  ^n-1  =  T  n 

■>  .  VT-*- 
-RT*N  G 
=  n 

where  the  matrix  operator  B  is 


B±  =  At2 (C  Q+RTN1 )V2 +1 

B.  SPECTRAL  FORMULATION 


(2.37) 


The  eguations  (2. 10) - (2. 15)  are  represented  spectrally 
in  the  horizontal.  The  variables  are  represented  as 
follows: 


m,  x,  o,  t)  =  J__ |m,C™(a.  t)p”(x)ei“A 


(2.  38) 


J  J 

Z  -Eii  C  Y 
m=-J  n=  m  n  n 


m  -m 


»  m  x  ill 

where  C  rs  some  variable  and  (Cn  )  =  (-1)  Cn 

m  =  zonal  wavenumber 

n  =  meridional  index,  and  n-Jml  gives  the  number  of  zeros 
between  the  poles  (-1  <  x  <  1)  of  the  associated  Legendre 
function 

J  =  truncation  limit  (for  this  study  triangular  truncation 
was  used  with  J=3) 

(1-1)  /2  nondimensional  zcnal  coordinate  index 
(1  <  1  <  16) 


I,  m.-  V  *  *.  *  *  ■».  ■  *,  -  •  \ •  • 

*  •  •  -  V'-  V.  t. 


A.*  w\.‘  A, 


/*  „*  A.'/*  .’v*’  -**  A  '  -  ■  A  ■>  A  A  A  A  . 


Note  that  the  separation  is  such  that  the  coefficients  c™ 

are  functions  of  time  and  the  vertical  and  the  spherical 
.  .  m  , 

harmonic  Y n are  horizontal  functions  of  space.  The  normali¬ 
zation  and  orthogonality  properties  of  the  Ym  allow  the 
.  n 

coefficients  to  be  obtained  as  follows: 


C™  =  <C,Ym> 
n  ’  n 


1  .  2  7T 
4tt  V 


f t1CYmdxdX 
- 1  n 


(2.39) 


The  non-linear  terms  are  computed  using  the  transform  method 
following  Haltiner  and  Williams  (1980).  The  longitudinal 
direction  is  treated  with  a  Fast  Fourier  Transform  and  the 
latitudinal  direction  uses  Gaussian  Quadrature.  The  number 
of  latitudes,  N  and  longitudes  U  satisfy. 

N  >  3J/2  +1  M  >  3J+1 

The  number  of  points  are  chosen  so  that  there  will  be  no 
aliasing  from  the  product  terms.  For  this  study  N=4  and 
M=16.  It  was  discovered  just  before  this  study  was 
completed  that  number  of  Gaussian  latitudes  used  to  compute 
the  latitudinal  integrals  was  actually  J*1  and  not  3J/2  +  1. 
However,  model  integrations  using  J  +1  latitudes  are  not 
qualitatively  different  than  these  using  3J/2  ♦  1  latitudes. 
Even  the  quantitative  differences  are  small,  and  will  occur 
mostly  in  the  non-linear  integrations.  Since  the  cases 
examined  in  this  study  are  carried  out  using  a  linear  form 
of  the  model,  one  would  not  expect  much  error. 
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III.  EXPERIMENT  DESIGN 


The  intent  of  this  study  is  to  specify  a  known  analytic 
heating  function  on  a  planetary  scale  and  then  integrate  the 
equations  until  the  model  atmosphere  reaches  a  steady  or 
guasi-steady  state.  These  steady  state  solutions  will  be 
altered  and/or  the  heating  will  be  changed.  The  equations 
will  then  be  integrated  again  to  determine  the  sensitivity 
of  planetary  waves  to  changes  in  the  forcing  and  initial 
conditions 

Heating  is  introduced  into  the  model  via  the  thermody¬ 
namic  eguation  (2.4).  Q  is  specified  to  be  a  Legendre 
function  of  the  form: 


Q(A,  x,  a)  =  i  E,  ,A™(a)P™(x)eimX 
m-  - n  -  |  m  J  n  n 

A™  is  specified  to  have  the  vertical  structure: 


(3.1) 


Am  = 


(3.2) 


Only  one  wave  form  is  forced  at  a  time  and  the  other  heating 
amplitudes  are  set  to  zero.  The  constant  A  is  chosen  to 
give  a  reasonable  temperature  response.  In  most  cases  a 
temperature  response  of  a  few  degrees  is  used. 

Before  the  heating  is  "turned  on"  the  model  atmosphere 
is  specified  to  be  in  a  state  of  solid  rotation 

U  =  cos 
V  =  0 

such  that  U  *  20  m/s(15.5  °/day)  at  the  equator.  The 
surface  pressure  field  is  chosen  to  balance  the  mean  flow. 
The  initial  vertical  temperature  distribution  is  given  by 
the  standard  atmosphere  values  consistent  with  the  pressure 


distribution  given  by  solid  rotation.  A  lapse  rate  of  6.5° 
per  kilometer  is  specified  tc  the  top  of  the  model 
atmosphere. 

TABLE  I 

Haves  Forced  for  linear  Cases 
Case 

1 

2 

3 

4 

5 

6 

In  all  of  the  experiments  the  model  is  run  in  a  linear 
mode.  In  this  mode  the  amplitudes  of  all  zonal  wavenumbers 
except  those  being  forced  are  held  constant.  The  model  is 
run  in  a  linear  mode  in  these  experiments  to  allow  a  more 
definitive  interpretation  of  the  results.  Six  linear  cases 
are  integrated,  each  case  with  a  different  forced  wave 
form.  For  example  in  case  1  (see  Table  I)  zonal  wavenumber 
one(m=1),  meridional  index  one  (n=1)  is  forced.  The  results 
from  these  cases  are  fairly  similar  so  only  the  longest  wave 
case  (case  1)  is  analyzed  in  detail. 

Two  classes  of  cases  are  examined  using  the  steady  state 
values.  In  the  first,  the  magnitude  of  the  heating  function 
is  changed  and  the  model  is  integrated  from  the  steady  state 
values.  In  the  second,  the  lagnitude  of  the  steady  state 
values  are  changed  and  then  used  as  initial  conditions  for 
further  model  integrations.  In  the  these  cases  the  heating 
function  is  the  same  as  the  one  originally  used  to  bring  the 
model  to  a  steady  state. 


Have  Forced 


m 

1 

1 

1 

2 

2 

3 


n 

1 

2 

3 

2 

3 

3 
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IV.  BESELTS 


Due  to  the  similarity  of  the  six  linear  cases  only  cases 
1,  5  and  6  will  be  discussed.  These  three  cases  represent 
the  three  zonal  wavenumbers ( 1 ,  2  and  3)  and  the  two  basic 
meridional  structures  possible  in  the  model.  Cases  1  and  6 
are  symmetric  in  their  meridional  structure  while  case  5  is 
anti-symmetric.  The  primary  tool  used  to  analyze  these 
cases  was  the  harmonic  dial.  These  diagrams  are  polar  plots 
of  phase  and  amplitude  for  a  given  spherical  harmonic.  This 
plot  has  the  advantage  that  the  guasi-sta tionary  and  tran¬ 
sient  nature  of  the  waves  are  easily  seen.  Fig.  4.1  is  a 
good  example  upon  which  the  merits  of  the  harmonic  dial  can 
be  explained.  In  this  figure  the  phase  and  amplitude  for 
the  temperature  wave  for  case  1  are  plotted  every  12  h.  The 
elapsed  time  since  the  start  of  the  integration  is  indicated 
every  5  days.  For  an  ideal  case  where  there  is  a  stationary 
and  a  transient  component  one  wculd  expect  to  see  a  circular 
pattern  which  is  confined  to  a  particular  quadrant.  This  is 
because  as  the  transient  wave  moves  through  the  stationary 
wave  it  will  come  in  and  out  cf  phase  with  the  stationary 
wave.  The  maximum  amplitude  occurs  when  the  transient  and 
stationary  waves  are  completely  in  phase  and  the  minimum 
amplitude  occurs  when  the  waves  are  180°  out  of  phase.  A 
rough  estimate  of  the  phase  speed  of  the  transient  wave  can 
be  obtained  by  observing  the  change  in  phase  with  time  as 
plotted  on  the  diagrams.  The  direction  of  propagation  can 
be  determined  by  observing  the  direction  of  the  phase  change 
on  the  dial.  A  counterclockwise  change  in  phase  indicates 
eastward  propagation,  a  clockwise  change  indicates  westward 
propagation.  In  Fig.  4.  1  the  phase  speed  of  the  temperature 
wave  at  a  =  .925  can  be  easily  estimated.  The  period  of  the 
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wave  is  21.5  days, 
is  360/21.5  =  16.74. 


thus  the  phase  speed  in  degrees  per  day 
The  direction  of  propagation  is  east- 


PHASE  VS  AMPLriCE 
FOR  TEMPERATURE 


Figure  4.1  Harmonic  Dial  for  Case  1  Temperature  Component. 

ward.  Another  important  point  to  note  here  is  that,  given 
the  uniformity  of  the  change  in  phase,  the  phase  speed  of 
the  transient  component  is  constant.  The  harmonic  dials  for 
surface  pressure,  temperature  and  the  zonal(U)  and  meri¬ 
dional  (V)  wind  components  for  the  lowest  model  level  are 
shown  in  Figs.  4.2  -  4.7  .  Only  those  components  which  have 
significant  amplitude  are  shown  in  these  figures.  In  all 
cases  the  pressure  and  winds  waves  show  a  much  more  irreg¬ 
ular  pattern  than  the  temperature  wave.  These  irregulari¬ 
ties  might  be  due  to  the  presence  of  a  number  of  transient 
modes  components  with  periods  much  less  than  that  indicated 
by  the  overall  pattern.  Table  II  gives  the  estimated  period 
(in  days)  and  phase  speed  {in  °/days)  for  the  temperature, 
surface  pressure  and  the  V  and  0  wind  component  waves  for 
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PHASE  VS  AMPLITUDE 
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Figure  4.2  Harmonic  Dials  for  Temperature 
and  Surface  Pressure  for  Case  1. 
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FUR  V  COMPONENT  OF  WIND  FOR  V  COMPONENT  OF  WIND 


Figure  4.3  Harmonic  Dials  for  Meridional  (0) 
and  Zonal  (0)  Hind  Case  1. 
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Figure  4.5  Harnonic  Dials  for  Meridional  (0) 
and  Zonal (0)  Bind  Case  5. 


Figure  4.6  Saae  as  Fig  4.4  for  Case  6 


PHASE  VS  AMPLITUDE 
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PHASE  VS  AMPLITUDE 
FDR  U  COMPONENT  OF  WIND 


the  three  cases  being  considered.  These  estimates  were 
obtained  from  Figs.  4.1  -  4.1  .  Phase  speeds  were  not 
calculated  for  the  pressure  and  wind  waves  for  case  5  due  to 
the  highly  irregular  nature  of  those  dials.  The  phase  speed 
of  the  temperature,  pressure  and  wind  waves  are  the  same  for 
wave  component  (3,3) .  Such  is  not  the  case  for  the  other 
cases,  though  the  phase  speed  are  similar.  The  computed 
phase  speed  of  the  different  variables  for  case  1  varies 
from  15.6  °/day  for  P  to  18.5  °/day  for  0.  The  other  cases 
have  similar  phase  speeds  for  the  different  components.  The 
general  pattern  which  emerges  from  the  data  is  that  the 
phase  speed  of  the  transient  components  seems  to  be  indepen¬ 
dent  of  wavenumber  and  that  the  waves  are  moving  at  approxi- 


TAB1E  II 

Period  (p) ,  and  Phase  Speed (c)  for  Transient  Components 


Case  Case  Case 

T~  ”5"  ~~5~ 

P  c  p  c  p  c 

P  22.  5  15.6  *  *  7f 5  16.7 

T  21.5  16.7  12.5  16.7  7.5  16.7 

U  19.5  18.4  *  *  7.5  16.7 

V  21.5  16.7  *  *  7.5  16.7 


mately  the  speed  of  the  mean  wind  (15.5  °/day)  .  To  test  this 
hypothesis  additional  integrations  were  accomplished  holding 
all  other  factors  the  same  but  increasing  the  speed  of  the 
mean  wind.  The  results  of  these  integrations  showed  an 
increase  in  phase  speed  of  the  component  waves  corresponding 
to  the  increase  in  the  mean  wind  speed 


A.  STEADI  STATE  CASE 


Since  the  purpose  of  this  study  was  to  see  how  sensitive 
the  steady  state  solutions  are  to  errors  in  the  forcing  and 
initial  conditions  it  is  important  that  the  steady  state 


solutions  of  the  model  be  examined  to  insure  that  they  are 
indeed  steady  and  that  these  solutions  are  consistent  with 
known  theory. 

To  obtain  steady  state  solutions  the  basic  equations 
were  integrated  up  to  30  days.  It  is  obvious  from  Figs. 
4. 1-4.7  that  after  a  rather  lengthy  integration  the  model 
atmosphere  has  yet  not  reached  a  steady  state  and  that  in 
fact  the  transient  modes  identified  in  the  last  section  are 
clearly  evident  throughout  the  period  of  integration. 
Though  there  are  indications  that  the  solutions  are  'lowly 
converging  as  indicated  by  the  begining  of  spiral  patterns 
in  Figs.  4. 1-4. 7  a  final  steady  state  would  not  be  reached 
for  some  time.  However,  by  averaging  over  the  period  of  the 
transient  wave,  it’s  effect  can  be  averaged  out.  The 
remaining  values  are  the  steady  state  values.  The  above 
procedure  was  performed  on  case  1.  To  insure  that  the  part 
of  the  solution  which  remained  after  averaging  was  truly  the 
steady  state  solution  another  integration  was  performed 
using  the  averaged  values  as  the  initial  conditions. 
Harmonic  dials  were  again  constructed  for  each  of  the  vari¬ 
able  wave  components.  If  the  fields  obtained  by  averaging 
were  truly  the  steady  state  solutions,  one  would  expect  to 
see  no  phase  or  amplitude  change  with  time,  i.e.  the 
harmonic  dial  would  collapse  to  a  single  point.  Figs.  4. 8 
and  4.9  contain  the  harmonic  dials  from  a  20-day  integration 
initialized  with  the  averaged  values  for  case  1.  The 
temperature  wave  is  almost  completely  stationary,  but  the 
pressure  and  velocity  components  have  small  growth  and  move¬ 
ment.  The  fact  that  the  pressure  and  wind  component  waves 
are  not  completely  stationary  is  another  indication  that 
there  are  transient  modes  in  these  waves  which  are  not 
contained  in  the  temperature  wave.  Consequently,  the  aver¬ 
aging  which  was  performed  over  the  period  of  the  temperature 
wave  did  not  average  out  all  of  the  transient  parts  of  the 
pressure  and  wind  waves. 
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It  appears  from  the  above  integrations  that  solutions 
obtained  by  averaging  over  the  period  of  the  transient 
component  are  very  close  to  steady  state  values.  In  order 
to  better  understand  these  solutions  they  will  be  compared 
to  the  steady  state  solutions  from  a  simplified  analytic 
model.  This  model  will  be  derived  in  the  next  section. 


Z  =  -p/p 


Eg.  4.1  is  the  continuity  equation.  Eg.  4.2  is  the  first  law 
of  thermodynamics  and  Eg.  4.3  is  the  geostrophic  relation. 
These  equations  were  obtained  ly  scaling  the  general  taro- 
clinic  equations  (Haltiner  and  Williams,  1980)  with: 
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1/a  ~1  R  <.01  e  -  100  fla/f  ~1 

Where  L  is  the  horizonal  scale,  R  is  the  Rossby  Number  and  e 
is  a  rotational  Froude  number.  Assuming  a  steady  state 
atmosphere  (|_  =  g  )  in  solid  rotation  (U=  wa  cosf'  )  and 
droping  ther’(Z)  term,  because  it  is  small.  Eg.  4.2  becomes: 


CD 


U'+ri  = 


3A  3 Z 


Q  ( Z ) 


(4.4) 


The  meridional  and  zonal  components  of  the  geostrophic  wind 
in  spherical  coordinates  are: 


y  _  _ 1 _  3  $ 1  (4.5) 

2aftcosfsint  3X 

U  -  1  H’  (4.6) 

2a2ft  sinY  3^ 

Using  Egs.  4.5  and  4.6  the  divergence  in  spherical  coordi¬ 
nates  can  be  written: 


_  2ftcosf 


8  V 
f 


2alsin4' 

Now  Eg.  4.1  can  be  written: 


11  -z-  ,  }  11’  -  n 

3 Z  u  a^ftsin2'*'  3X  "  u 


Assuming  solutions  of  the  form: 


=  <KZ)elmAP™('0 


Z  =  W( Z)elmApm( f) 
n 


Q  =  N(Z)eimAP™(Y) 


Eg.  4.4  and  4.8  become: 


cuim-M  Z)  +TW(  Z)  =  N(Z) 
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(4.7) 


(4.8) 


(4.9) 
(4.  10) 
(4.  11) 


(4.  12) 


W'-W- 


im* ,  .  0  (4.  13) 

a  2  Q,  s  i  n 2  'F 

Eg.  4.12  and  4.13  can  be  combined  to  obtain  a  single  equa¬ 
tion  for  W 


oja2ft  sin2 '¥  ( W" -W '  )  +  TW  =  N  (4.14) 

Eg.  4.14  is  just  a  2nd  order  n cnhomogenous  ordinary  differ¬ 
ential  equation  for  W  whose  general  solution  is  of  the  form: 


W  =  C.w.+C-W.+W 

11  12  p 


Where  •  1*2  are  the  solutions  to  the  homogeneous  part  of 

Eg.  4.14  and  W  p  is  the  particular  solution.  Assuming  N  (Z)  = 


N0e-Ct  7  the  solution  to  Eg.  4.14  is 


W  =  C1er*“*C2er-ZtAe_a  Z 


=  eZ/2(C1e'l'Z*c2e-<3Iz, 


where 

=  j  -»,'l/4-r/oja2r2sin'l'  =  i  +q2 

«2=  seme  positive  number 


(4.  15) 


A  =  NQ/ (coa2  sin2  TC  a1* +a2  ) +T  ) 

There  are  two  possible  solution  cases  for  the  homogeneous 
part  of  Eg.  4.14.  If  y2  is  positive  the  solution  consists 
of  an  exponentially  growing  and  decaying  part.  After 
applying  the  boundary  conditions  w[0)=0  and  finite  energy  at 
Z=®#  i.e.  (  pHI  2  )  =  e-7!^)2  the  solution  becomes 


W(Z)  =  A(e 

An  exfression 
4. 16  into  Eg. 


Z  Z/2  -q 2  Z . 

-e  e  n  ) 

for  t  (z)  can 

4.  13 


(4.  16) 

be  obtained  by  substituting  Eg. 
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<KZ) 


(4.  17) 


=  -a2«sin24'AterZ(r-l)+e'a2Z(l+  a2)} 

If  q2  is  negative  the  solution  will  be  oscillatory  in  Z  and 
of  the  form 

W(Z)  =  eZ/2(C1e"iyZ+C2feluZ)+Ae"a2z  l4-18) 

where 

U  =  /r/ti)a2^s"in2f-  1/4  (4-  19) 

After  applying  the  radiation  boundary  condition  at  Z=“i.e. 
w(Z)  <t  (Z)  >0  and  the  lower  boundary  condition  w  =  0  at  Z  =  0  Eg. 
4.  18  becomes 

WCZ)  =  A{(|i-M)e(1'5+ilj)Z-i(-l-a2)e-a2Z}  t4*20) 

The  corresponding  expression  for  (z)  is 

$  (  Z )  =  a2fisin2 'I'AC  e(  15  +  111 )  Z-i( -l-a2  )e-a  Z}  f4--^1) 

Egs.  4.16,4.17  and  Egs.  4.20,4.21  can  be  substituted  into 
Eqs.  4.9  and  4.10  to  obtain  expressions  for  <f>and  Z.  After 
taking  the  real  part  the  g2  >  0  solution  are 

<P'  =  a2ftsin2YA{erZ(r-l)+e~a2Z(l  +  a2  )}  sin(mA  1  Pm 

n 

y  . ,  -a2Z  rZ ,  ,  , . m 

Z  -  A(e  -e  )cos(mA)P 

n 

The  q2  <  0  solutions  are 

<t> '  =  a2^  s in2  f  A{  eZZ  2  ( ysinC  yZ+mA  )  (4.24) 

+  MCOS  (  u  Z  +  mA  )  +  (  1  +  a2  )  e  a  ^sinCmA)}?111 

n 


(4. 22) 

(4.23) 
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Z  =  A(  e 


cos(pZ)  )cos(m.\) 


(4. 25) 


-a2Z  ,  Z/2 
-  ( e 


+  (  e~'  2sin(  ijZ)  )sin(mA)  pm 

If  q2  >  0  the  response  of  the  atmosphere  to  heating  is 
trapped  near  the  surface  and  the  vertical  structure  is  one 
in  which  the  waves  have  no  tilt  with  height.  If  g2  <  0  the 
perturbation  due  to  the  heating  will  radiate  to  great 
heights  and  the  waves  will  tilt  westward.  In  both  cases  the 
steady  state  represents  a  balance  between  advection,  adia¬ 
batic  warming  (coolin g)  and  diatatic  heating  (cooling)  .  The 
longitudinal  and  latitudinal  structure,  for  each  case,  will 
be  determined  by  the  forcing  specified.  That  is,  if  waven¬ 
umber  one  is  specified  in  the  forcing  term,  the  steady  state 
solution  will  have  a  wave  one  form-  In  addition,  since  the 
winds  are  geostrophic  their  structure  will  be  determined  by 
the  structure  of  <6*.  Eq.  4.24  can  be  written 


4> '  =  -A{  ( ^eZZ 2  sin(  u  Z) -eZZ  2y  cos  (  y  Z )  )cos(mA) 

+(feZ/2cos(yZ)-yeZ/2sin(uZ)-(l+a2 ) e~a2 Z ) sin ( m) ) } Pm 

L  n 

which  can  be  written  in  the  form 


{ G( Z ) cos ( nA  +H ( Z )  ) }  A(  )  Pm(  V )  1  J 

n 

This  form  allows  comparison  with  the  phase  and  amplitude  of 
the  model  Legendre  functions.  If  the  model  steady  state 
solutions  really  are  steady  state  one  would  expect  them  to 
at  least  qualitatively  agree  with  the  analytic  solutions 
derived  above.  The  differences  between  the  numerical  model 
atmosphere  and  the  analytic  atmosphere  are: 

1.  A  nondimension  vertical  coordinate  Z  is  used  in  the 
analytic  model  while  sigma  coordinates  are  used  in 
numerical  model, 

2.  Friction  is  not  included  in  the  analytic  model,  and 
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3. 


r(Z)  is  assumed  to  be  constant  in  the  analytic  model 
while  it  is  not  in  the  numerical  model. 

T(Z)  is  not  constant  in  the  numerical  model  because  r(z)  is 
constant.  The  fact  that  T{Z)  is  not  constant  given  r  (z)  is 
can  be  shown  by  examining  the  eguation  for  r  (Z) 


r(z) 


a  + 

Cp 


1  3T 
H  9  Z 


(4. 28) 


After  converting  to  z  coordinates  Eg.  4.28  becomes 


H  (fi  + 

T  ^Cp 


3T, 

3Z' 


(4. 29) 


Since  the  lapse  rate  is  a  constant  6.5  °/Km  in  the  model  the 
second  term  in  Eg.  4.29  is  constant  but  T  decreases  with  z 
so  the  first  term  is  not  constant  thus  r  (Z)  will  not  be 
constant  in  the  model  atmosphere.  However,  since  the 
decrease  in  temperature  is  rather  small  compared  to  the 
mean  temperature  T(Z)  is  nearly  constant  so  that  qualita¬ 
tive  comparisons  are  still  possible.  In  addition,  the 

vertical  structure  of  the  forcing  in  the  analytic  model  is  a 

Z 

simple  exponential  form  e  while  in  the  model  it  is  of 

the  form  ctt  =  e?^s  .  However,  given  that  is  chosen  such 

300  -a2Z 

that  the  e-folding  depth  of  e  is  the  same  as  that  for 
P-Pc 

e300  t^ien  two  functions  will  be  similar  except  for  the 

upper  levels  where  both  functions  are  small  and  have  little 
effect  on  the  total  solution.  Fig.  4.10a  is  a  plot  of  the 
numerical  model  forcing  vertical  structure  where  Pg  is 
assumed  to  be  constant  and  egual  to  lOOOmb.  This  is  not  a 
bad  assumption  because  the  model  contains  no  topography  and 
the  maximum  pressure  perturbations  are  only  a  few.  millibars. 
This  assumption  also  allows  gualitative  comparisons  between 
variables  in  and  P  coordinates.  The  e-folding  depth  here 
occurs  at  P=700mb.  For  P=700  ,  Z=-ln  (700/1000)  =-ln  (.  7)  . 

Thus,  to  have  the  sane  folding  depth  a*  must  be  choosen  to 
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be  equal  to  1/(-ln.7).  Fig.  4.10b  is  a  plot  of 
converted  to  P  coordinates.  Note  that  the  plots  are  very 

similar  up  to  400mb.  Figs.  4.11a  and  4.11b  are  plots  of 
G(Z)  and  H  (Z)  vs  pressure  using  the  following  estimate  for 
the  various  parameters  in  Eg.  4.19:  w  =(20  m/s)/a,  r~  (Z) 
H2N2,  H=10km,  N2  =  1 0~* s~ 1 ,  a  =  6.  37  x  10‘  m  and  &  =  7.292  x 

10-ss-»  .  Figs.  4.12  and  4.13  contain  plots  of  the  steady 
state  model  solutions  for  the  meridional  (V)  and  zonal  (U) 
wind  vs  sigma  for  case  1.  Since  the  winds  are  geostrophic 
the  vertical  structure  U  and  V  can  be  compared  with  the 
vertical  structure  of  $.  Note  the  striking  similarities 
between  the  two  sets  of  figures.  Both  amplitude  plots  show 
a  slight  decrease  in  the  lower  layers  and  an  exponential 
increase  near  the  top  of  the  atmosphere.  Also,  given  Figs. 
4.14,  4.15  and  4.16,  it  is  obvious  that  the  model’s  steady 
state  solutions  have  the  longitudinal  wave  one  form  speci¬ 
fied  in  the  case  1  forcing  term.  These  figures  are  hori¬ 
zontal  cross-section  for  vorticity,  meridional  velocity  and 
temperature  near  the  top  and  bottom  of  the  model  atmosphere. 
Contour  intervals  are  indicated  at  the  top  of  each  cross- 
section.  The  vertical  scale  at  the  bottom  of  the  figures  is 
used  to  represent  topography,  which  is  not  included  in  this 
model.  It  is  clear  that  the  steady  model  solutions  are 
qualitatively  consistent  with  the  analytic  solutions  and 
that  the  method  used  to  determine  these  values  was  a  sound 
one.  The  question  might  now  be  asked  why  are  the  model 
solutions  oscillatory  in  nature?  Is  it  not  possible  that 
they  could  be  of  the  exponential  decaying  type?  These  ques¬ 
tions  can  be  answered  by  examining  Eg.  4.19  more  closely. 
Given  that 

r/a>a2sin24'  >  ~  (‘*•^0) 
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Figure  4. 12  Model  Vertical  Amplitude 
and  Phase  Structure  for  Case  1. 
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the  solution  will  be  oscillatory.  In  the  model,  <o  is  speci¬ 
fied  so  that  the  zonal  wind  is  equal  to  20  m/s  at  the 
equator  so  that  U=  gj  a  cos  =  U=<ua  at  =0  thus  o)  =  U/a.  Given 
T~H2  N2,  H ~1 0  km,  N2  =10-*s~i,  a=6.37  x  10*  m,  Q  =7.292  x 

10-5s_l  one  can  try  to  solve  f cr  a  latitude  where  the  solu¬ 
tion  will  no  longer  be  oscillatory.  After  using  the  above 
values  Eg.  4.30  becomes 


sinf  =  20.84 


(4.31) 


which  has  no  solution.  Thus  there  is  no  point  where  the 
solution  would  not  be  oscillatory.  To  obtain  non- 

oscillatory  solutions  one  would  have  to  increase  the  zonal 
wind  speed  to  almost  100  m/s.  Using  these  steady  state 
solutions  it  is  now  possible  to  examine  the  sensitivity  of 
these  solutions  to  changes  in  the  forcing  and  initial 
conditions. 


C.  SENSITIVITY  EXPEBIHENTS 

Using  the  steady  state  solutions  from  case  1  as  initial 
conditions  two  additional  experiments  are  carried  out  to 
determine  the  sensitivity  of  the  planetary  waves  to  'errors' 
in  the  forcing  or  initial  conditions.  In  the  first  experi¬ 
ment  errors  are  introduced  in  the  forcing  field.  In  fact  a 
100%  error  was  introduced  by  turning  off  the  forcing  alto¬ 
gether.  The  model  is  then  integrated  out  to  20  days.  Figs. 
4.17  and  4.18  contain  the  harmonic  dials  for  the  above  inte¬ 
gration.  Note  that  for  the  most  part  the  once  quasi¬ 
stationary  solutions  are  now  almost  purely  transient  since 
the  phase  through  the  entire  range  of  phase.  In  the  previ¬ 
ously  observed  pattern  the  stationary  plus  transient  nature 
of  the  waves  was  indicated  by  a  change  of  phase  and  ampli¬ 
tude  which  occured  in  a  particular  quadrant  of  the  dial. 


The  estimated  phase  speeds  obtained  from  Figs  4.17  and  4.18 
show  these  transient  waves  to  again  be  moving  at  the  speed 
of  the  mean  wind.  The  fact  that  the  waves  are  pure  tran¬ 
sient  now  is  not  surprising  since  the  heating  is  no  longer 
present  to  fix  the  wave  to  a  particular  location. 

In  the  second  experiment  errors  are  introduced  into  the 
rotational  part  of  the  steady  state  wind  by  reducing  those 
components  by  20 %.  The  model  is  again  integrated  out  to  20 
days.  Figs.  4.19  and  4.20  ccntain  the  harmonic  dials  for 
this  integration.  From  these  figures  one  can  see  that  the 
only  wave  which  is  moving  eastward  is  the  temperature  wave. 
All  the  other  component  waves  are  moving  westward.  The 
temperature  wave  has  a  phase  speed  of  about  12  °/day  at  the 
equator  which  is  about  25%  slower  than  the  mean  wind.  It  is 
evident  that  the  j  are  a  number  of  high  frequency  components 
present  in  the  other  component  waves  which  are  not  present 
in  the  temperature  wave.  Since  there  is  no  simple  circular 
pattern  for  the  wind  and  pressure  waves  it  is  difficult  to 
obtain  the  period  of  the  transient  part  of  the  solution. 
The  exact  phase  speed  of  these  wave  is  not  important.  The 
important  thing  to  note  here  is  that  the  waves  are  moving 
slowly  westward. 
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PHASE  VS  AMPLITUDE 
FDR  SURFACE  PRESSURE 


4.17  Same  as  Fig  4.2  fcr  the  Case  with  No  Heating 
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fhase  vs  amplitude  phase  vs  amplitude 

FDR  TEMPERATURE  FOR  TEMPERATURE 


Figure  4.19  Same  as  Jig  4.2  for  the  Case 
with  Botational  Hind  Errors. 


V.  PI SC0S5I0H  AHE  CONCLUSIONS 


The  results  of  the  previous  section  are  best  discussed 
in  terms  of  normal  mode  theory.  A  number  of  investigators 
(Geisler  and  Dickinson,  1976,  Daley  et  al,  1981,  and 
Kasahara  1976)  have  used  this  method  to  examine  the  tran¬ 
sient  part  of  planetary  waves.  The  method  basically 
consists  of  solving  the  linearized  set  of  eguations 
describing  oscillations  of  a  stratified  resting  atmosphere 
on  a  spherical  earth  by  separating  the  eguations  into  a 
latitudinal  structure  equation  and  vertical  structure  equa¬ 
tion.  The  vertical  structure  equation  determines  the 
vertical  modes  of  oscillation.  The  separation 
constants (eigenvalues)  from  the  vertical  structure  equation 
can  be  used  to  find  the  eigenfunctions  of  the  latitudinal 
structure  equation.  The  latitudinal  structure  equation  for 
a  given  eigenvalue  is  identical  to  a  fundamental  equation  of 
a  free  surface,  if  the  eigenvalue  is  replaced  by  the  uniform 
depth  of  the  fluid.  This  equation  was  first  obtained  by 
Laplace  in  his  study  of  free  oscillations  of  shallow  water 
over  a  rotating  sphere.  For  this  reason  the  horizontal 
structure  equation  is  often  referred  to  as  Laplace's  tidal 
equation  and  the  eigenvalues  are  referred  to  as  equivalent 
depths.  The  mode  associated  with  the  largest  equivalent 
depth  is  referred  to  as  the  external  or  barotropic  mode. 
This  is  because  this  mode  has  a  vertical  profile  of  hori¬ 
zontal  divergence  with  the  same  sign  throughout  the  atmos¬ 
phere.  The  other  modes  associated  with  smaller  equivalent 
depths  are  referred  to  as  internal  modes.  These  modes  are 
oscillatory  in  the  vertical.  For  a  given  equivalent  depth, 
zonal  wavenumber  and  meridional  mode  there  exist  three  sets 
of  eigenfrequencies  and  eigenfunctions.  Two  are  eastward- 


and  westward-propagating  gravity  waves  and  the  third  is  a 
westward  propagating  Rossby-Haurwitz  type  wave.  The  phase 
speeds  of  these  waves  are  related  to  eguivalent  depth,  meri¬ 
dional  mode  and  zonal  wavenumber.  For  the  Rossby  type  wave 
the  fastest  phase  s  peed  occurs  with  the  largest  eguivalent 
depth  and  smallest  meridional  mode  and  zonal  wavenumber. 
The  period  of  this  gravest,  symmetric  zonal  wavenumber  one 
Rossby  mode  is  almost  exactly  5  days.  There  is  observa¬ 
tional  evidence  for  the  existence  of  this  wave (Madden  and 
Julian,  1972)  and  there  is  also  evidence  that  some  of  the 
planetary  wave  error  in  numerical  models  is  in  part  due  to 
the  spurious  excitation  of  these  large  scale  external  Rossby 
modes  (Daley  et  al,  1981).  The  following  transient  modes 
which  were  observed  in  the  three  different  parts  of  this 
study  are: 

1.  the  mode  which  moved  at  the  speed  of  the  mean  wind 

in  the  case  where  the  model  atmosphere  was  being 
brought  to  a  steady  state, 

2.  the  mode  which  moved  at  the  speed  of  the  mean  wind  in 
the  case  where  the  model  heating  was  turned  off, 

3.  the  westward  moving  modes  in  the  case  where  error  was 
introduced  into  the  rotational  part  of  the  wind. 

These  modes  may  be  explainable  in  terms  of  these  external 
and  internal  Rossby  modes.  The  mode  which  is  moving  at  the 
speed  of  mean  wind  may  be  a  slow  westward  moving  internal 
mode  which  is  overwhelmed  by  the  mean  wind.  The  westward 
moving  modes  in  the  last  experiment  are  most  likely  internal 
modes  as  well  but  are  probably  associated  with  larger  eguiv¬ 
alent  depths  than  those  of  the  earlier  experiments.  Since 
phase  speed  for  a  given  wavenumber  and  meridional  mode 
increases  with  increasing  eguivalent  depth,  these  waves 
would  have  a  greater  westward  phase  speed  and  would  be  less 
counterbalanced  by  the  mean  wind.  Given  a  large  enough 
westward  phase  speed  one  would  actually  see  a  westward 
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propagation  indicated  on  the  harmonic  dials.  The  low 
frequency  modes  in  the  first  two  parts  of  this  study  were 
most  likely  excited  by  the  discontinuous  nature  of  the 
forcing.  Numerical  studies  by  Geisler  and  Dickinson  (1976) 
indicate  that  even  when  the  forcing  is  gradually  turned  on 
that  low  frequency  Eossby  modes  are  excited.  An  analytic 
study  by  Clark  (1972)  also  showed  that  a  mode  which  moved  at 
the  speed  of  the  mean  wind  would  be  excited  by  a  switch  on 
of  a  vertical  velocity  or  temperature  disturbance.  However, 
this  mode  was  only  cne  of  a  number  of  modes  which  would  be 
excited  and  it  was  also  found  to  decay  with  time.  Thus,  it 
is  possible  that  the  transient  modes  which  were  observed  to 
move  at  the  speed  of  the  meat  wind  are  in  fact  internal 
Eossby  modes.  The  fact  that  internal  Eossby  modes  might 
have  been  excited  when  errors  were  introduced  into  the  rota¬ 
tional  part  of  the  wind  is  consistent  with  the  results  of 
Daley  et  al  (1981)  who  found  that  errors  in  the  rotational 
wind  will  cause  a  spurious  excitation  of  these  modes. 
However,  it  must  be  pointed  out  that  since  the  results  of 
this  study  were  analyzed  using  spherical  harmonics  rather 
than  normal  modes  it  is  impossible  to  confirm  the  above 
hypothesis.  The  question  that  might  now  be  asked  is  why 
wasn't  the  external  barotropic  node  excited  in  these  experi¬ 
ments?  The  studies  by  Geisler  and  Dickinson  and  Clark  found 
that  this  mode  would  be  excited.  Clark  found  that  this  mode 
would  in  fact  be  the  dominate  mcde.  The  crucial  difference 
between  the  analytic  model  of  Clark  ,  the  numerical  model  of 
Geisler  and  Dickinson  and  the  model  in  this  study  is  that 
this  model  does  not  have  a  vertical  wind  shear  while  the 
aforementioned  models  do.  This  lack  of  vertical  shear  is 
the  most  probable  reason  why  the  barotropic  mode  was  not 
excited.  In  the  absence  of  a  mean  vertical  shear  the 
heating  specified  in  this  model  will  produce  low  pressure  at 
the  surface  and  a  high  pressure  aloft.  The  vorticities 


associated  with  each  will  be  of  the  same  magnitude  so  that 
the  vertically  averaged  vorticity  vanishes  and  there  can  be 
no  barotropic  mode  response.  Khen  a  mean  vertical  shear  is 
present  there  is  a  differential  advection  affect  and  the 
upper  and  lower  vorticities  dc  not  exactly  cancel,  thus 
allowing  the  possiblity  of  a  barotropic  mode 

Given  the  slow  phase  speed  of  the  transient  modes 
excited  in  this  study  by  ’errors'  in  the  forcing  and  initial 
conditions  it  does  not  appear  that  planetary  waves  are  very 
sensitive  to  errors.  However,  this  apparent  lack  of  sensi¬ 
tivity  is  most  likely  due  to  the  simplified  conditions  used 
in  the  model  rather  than  to  the  inherent  nature  of  planetary 
waves.  Specifically,  the  lack  of  wind  shear  prohibits  the 
excitation  of  the  fast  external  Rossby  wave.  The  external 
mode  has  a  fast  phase  speed  (  72  °/day)  and  if  excited  it 
could  rapidly  propagate  throughout  the  region  of  integra¬ 
tion. 

The  results  obtained  by  heating  in  this  study  might  also 
have  been  obtained  by  including  topography  along  with  a 
stronger  surface  friction  term.  Future  studies  should 
investigate  the  sensitivity  of  planetary  wave  to  errors  in 
topography  as  well  as  heating.  In  addition,  future  studies 
will  have  to  use  a  more  realistic  vertical  wind  profile. 
Due  to  the  lack  of  a  realistic  wind  profile,  which  prohibits 
the  excitations  of  the  external  barotropic  mode,  any  conclu¬ 
sions  about  the  sensitivity  of  planetary  waves  to  initial 
conditions  and  forcing  must  be  delayed  until  more  comprehen¬ 
sive  studies  are  completed. 


LIST  OF  REFERENCES 


Arakawa,  A. ,  and  M.  J.  Suarez,  1983:  Vertical  differencing 
of  the  primitive  equations  in  sigma  coordinates.  Mon. 
Sea.  Rev.,  Ill,  34-45.  - 

Baumhefner,  p. ,  and  P.  Downey,  1978:  A  comparison  of  six 

wintertime  forecasts  from  several  numerical  weather 
prediction  models.  Proc.  12th  Stanstead  Seminar. 
McGill  Oniversity  Publicaf  Ions"n-Meteorology  No.  " 12 1, 
98pp. (Department  of  Meteorology  McGill  Oniversity,  805 
Sherbrooke  St.  W.,  Montreal  P.Q.,  H3A  2K6  Canada.) 

Bourke,  W.,  1974:  A  multi-level  spectral  model  1. 

Formulation  and  hemispheric  integrations.  Mon.  Sea. 

Bey.,  102,  687-701.  -  - 

Clark,  J.  H.  E.,  1972:  The  vertical  progaga tion  of  forced 
atmospheric  planetary  waves.  J.  Atmos.  Sci.,  29, 
1430-1451  .  - -  - 

Daley.  R.,  J.  Tribbia,  and  D.  Rilliamson,  1981:  The  excita¬ 
tion  of  large-scale  free  Rossby  waves  in  numerical 
weather  prediction.  Mon.  Sea.  Rev. .  109,  183  6-1861. 

Geisler,  J.  E.  and  R.  E.  Dicksison,  1976:  The  five-day  wave 
on  a  sphere  with  realistic  zonal  winds.  J.  Atmos. 

Sci.,  33,  632-641.  - - 

Haltiner,  Gt  J.  and  R.  T.  Williams,  1980:  Numerical 

Prediction  and  Dynamic  Meteorology.  John  WileT- E“S6ng7 
TT77pp.  ~ 

Hoskins,  B.  J.  and  A.  J.  Simmons,  1975:  A  multi-layer  spec¬ 
tral  model  and  the  semi-implicit  method.  Ouart.  J. 
Boy.  Meteor .  So c . ,  101,  63  7-655. 

Kasahara,  A.,  1976:  Normal  modes  of  ultralong  waves  in  the 

atmosphere.  Mon.  Sea.  Rev.,  104,  669-690. 

Lambert,  S.,  and.P.  Merilees,  1978;  A  study  of  planetary 
wave  errors  in  a  spectral  numerical  weather  prediction 
model.  Atmos.  Ocean.  16,  197-211. 

Lorenz,  E.,  1969:  The  predictability  of  a  flow  which  poss¬ 

esses  many  scales  or  motion.  Tellus.  21,  248-307. 

Madden,  R.  A.  and  P.  A.  Julian,  1972:  Further  evidence  of 
qlotal-scale  5-day  pressure  waves.  J.  Atmos.  Sci.,  29, 
1464-1469.  - - - 


Morse,  P.  J.,  1983:  NOGAPS  verification  using  spectral 

components.  Masters  Thesis,  Naval  Postgraduate  School, 
71pp.  “  * 

Rosmond,  T.  E. ,  1977:  Personal  correspondence. 


INITIAL  DISTRIEOTION  LIST 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314 

2.  Library,  Code  0  142 
Naval  Postgraduate  School 
Monterey,  CA  93943 

3.  Commander 

Naval  Oceanography  Command 

N STL  Station 

Bay  St  Louis,  MS  39522 

3.  Commander 

Air  Weather  Service 

Scott  Air  Force  Base,  IL  62225 

4.  Commanding  Officer 

Fleet  Numerical  Oceanography  Center 
Monterey,  CA  93943 

5.  Commanding  Officer 

Air  Force  Global  Weather  Central 
Offutt  Air  Force  Base,  NE  68113 

6.  Cfficer-in-Charge 

Naval  Environmental  Prediction  Research 
Facility 

Monterey,  CA  93943 

7.  Prof.  C.  N.  K.  Mooers.  Code  63Mr 
Naval  Postgraduate  School 
Monterey,  CA  93943 

3.  Prof.  R.  J.  Renard,  Code  63Rd 
Naval  Postgraduate  School 
Monterey,  CA  93943 

9.  Major  Patrick  Herod 
AFJT/CIRF 

Wrlght-Patterson  AFB,  OH  45433 

10.  Air  Weather  Service  Technical  Library 
Scott  AFB,  IL  62225 

11.  Captain  Michael  D.  McAtee,  Code  63 
Naval  Postgraduate  School 
Monterey,  CA  93943 

12.  Director  of  Research  and  Administration 
Code  012 

Naval  Postgraduate  School 
Monterey,  CA  93943 

13.  Prof.  R.  T.  Williams,  Code  63Wu 
Naval  Postgraduate  School 
Monterey  CA  93943 


No.  Copies 
2 


2 


1 


1 


1 


1 


1 


1 


1 


2 


1 

2 


1 


3 


60 


14.  Prof.  M.  A.  Rennick,  Code  63 
Naval  Postgraduate  School 
Monterey  CA  93943 

15.  Dr.  T  Rosmond 

Naval .Environmental  Prediction  Research 
Facility 

Monterey,  CA  93943 

16.  Dr.  J.  Hovermale 

Naval . Environmental  Prediction  Research 
Facility 

Monterey,  CA  93943 

17.  Prof.  C.  H.  Wash,  Code  63 
Naval  Postgraduate  School 
Monterey,  CA  93943 


fvV. 


s>V- 

v'.-. 

v\-. 


k«*. 


